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We report first results on elliptic flow of identified particles at mid-rapidity in Au+Au colli- 
sions at ,/s;^ = 130 GeV using the STAR TPC at RHIC. The elhptic fiow as a function of 
transverse momentum and centrality differs significantly for particles of different masses. This 
dependence can be accounted for in hydrodynamic models, indicating that the system created 
shows a behavior consistent with collective hydrodynamical fiow. The fit to the data with a 
simple model gives information on the temperature and fiow velocities at freeze-out. 
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The goal of the ultra-relativistic nuclear collision pro- 
gram is the creation of a system of deconfined quarks 
and gluons pi . The azimuthal anisotropy of the trans- 
verse momentum distribution for non-central collisions is 
thought to be sensitive to the early evolution of the sys- 
tem H . The second Fourier coefficient of this anisotropy, 
V2, is called elliptic flow. It is an important observable 
since it is sensitive to the rescattering of the constituents 
in the created hot and dense matter. This rescattering 
converts the initial spatial anisotropy, due to the almond 
shape of the overlap region of non-central collisions, into 
momentum anisotropy. The spatial anisotropy is largest 
early in the evolution of the collision, but as the system 
expands and becomes more spherical, this driving force 
quenches itself. Therefore, the magnitude of the observed 
elliptic flow should reflect the extent of the rescattering 
at relatively early time [g| . 

Elliptic flow in ultra-relativistic nuclear collisions was 
first discussed in Ref. H and has been studied intensively 



in recent years at the Alternating Gradient Synchrotron 
(AGS) |, |, Super Proton Synchrotron (SPS) |, 0, i 
and at the Relativistic Heavy Ion Collider (RHIC) |£] 
energies. The studies at top AGS energy and SPS have 
found that elliptic flow at these energies is in the plane 
defined by the beam direction and the impact parame- 
ter, V2 > 0, as expected from most models. The pion 
elliptic fiow for relatively peripheral collisions increases 



with beam energy |10 from about 0.02 at the top AGS 
energy |h 0.035 at the SPS [| to about 0.06 at RHIC 
energies^ . The increased magnitude of the integrated 
elliptic flow at RHIC reaches the values predicted by hy- 
drodynamical models, which are based on the assumption 
of complete local thermalization. 

The first elliptic flow results from RHIC were for 
charged particles. The differential charged particle flow, 
V2{pt), shows an almost linear rise with transverse mo- 
mentum, pt, up to 1.5 GeV/c. At pt > 1.5 GeV/c, the 
i^2{pt) values start to saturate, which might indicate the 



onset of hard processes [|[ ^ |T^, |Tj. The behavior 
of V2{pt) up to 1.5 GeV/c is consistent with a hydrody- 
naniic picture. Hydrodynamics assumes complete local 
thermalization at the formation of the system, followed 
by an evolution governed by an Equation Of State (EOS) . 
However, the pt-integrated elliptic flow, V2, as a function 
of centrality, and the differential V2{pt), show little sen- 
sitivity to the EOS used |lj|. Studies of the mass depen- 
dences of elliptic flow for particles with pt < 1.5 GeV/c 
provide important additional tests of the hydrodynam- 
ical model |jl^. Similar to the identified single particle 
spectra, where the transverse flow velocity has been ex- 
tracted from the mass dependence of the slope parame- 
ter ijl^l , the V2 (pt ) for different mass particles allows the 
extraction of the elliptic component of the flow veloc- 
ity ||l^, Q. Moreover, the details of the dependence of 
elliptic flow on particle mass and transverse momentum 
are sensitive to the temperature, transverse flow veloc- 
ity, its azimuthal variation, and source deformation at 
freeze-out. In this Letter we report the first results for 
identified particle V2{m,pt) in Au-|-Au collisions at RHIC 
at y/s;^ = 130 GeV. 

The Solenoidal Tracker At RHIC (STAR) ^, is ide- 
ally suited for measuring elliptic flow due to its azimuthal 
symmetry and large coverage. The detector consists of 
several sub-systems in a large solenoidal magnet. The 
Time Projection Chamber (TPC) covers the pseudora- 
pidity range \r]\ < 1.8 for collisions in the center of the 
TPC. The magnet was operated at a 0.25 Tesla field, 
allowing tracking of particles with pt > 75 MeV/c. Two 
Zero Degree Calorimeters p^] located at | sm6\ < 0.002, 
which mainly detect fragmentation neutrons, are used in 
coincidence for the trigger. The TPC is surrounded by 
a scintillator barrel which measures the charged particle 
multiplicity within jryj < 1, for triggering purposes. 

For this analysis, 120,000 minimum-bias events were 
selected with a primary vertex position within 75 cm lon- 
gitudinally of the TPC center and within 1 cm radially of 
the beam line. For determination of the event plane [^ , 
charged particle tracks were selected with 0.1 < pt < 2.0 
GeV/c. All tracks used in this analysis passed within 
2 cm of the primary vertex and had at least 15 measured 
space points. Also, the ratio of the number of measured 
space points to the expected maximum number of space 
points for that particular track was required to be greater 
than 0.52, suppressing split tracks from being counted 
more than once. The tracks used for the determination 
of the reaction plane were within \ri\ < 1.0, and the tracks 
used to calculate the elliptic flow were within \ri\ < 1.3. 
These cuts are similar to the ones used in Ref. H, and 
the analysis results presented here are not sensitive to 
those cuts. 

The pions, protons, and antiprotons were selected ac- 
cording to their specific energy loss [dE/dx) in the TPC 
in the transverse momentum range of 0.175 - 0.75 GeV/c, 
0.5 - 0.9 GeV/c and 0.3 - 0.9 GeV/c, respectively. The 
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FIG. 1: Charged particle identification for central events at 
\rj\ < 0.1 for the momentum range of 380 - 400 MeV/c and 
700 - 720 MeV/c. The shaded areas show the selected ranges 
for the different particles. 



lower Pt values were chosen such that the energy loss in 
the detector was negligible. At low momentum, the pro- 
ton background due to secondary interactions with the 
detector material is significant. Therefore, only protons 
above a transverse momentum of 0.5 GeV/c were used 
in this analysis. At a momentum (p) of 0.5 GeV/c, the 
dE/dx resolution was found to be of the order of 11% 
for a typical long track in the STAR TPC. The kaons, 
because of their overlap in dE/dx with the electrons and 
positrons, were selected in the momentum ranges of 0.30 
- 0.40 GeV/c and 0.60 - 0.70 GeV/c. The raw yields of 
the pions, kaons, protons, and antiprotons were obtained 
from fitting the dE/dx distributions for each 77, p and cen- 
trality bin with multiple Gaussians, and requiring greater 
than 90% purity (particle's yield divided by the sum of 
yields for all particles at that same dE/dx). Figure |l| 
shows the dE/dx distributions for negative tracks in the 
TPC for central events and \rj\ < 0.1. The upper panel 
shows the dE/dx distribution for the momentum range 
of 0.38 - 0.40 GeV/c and the multiple Gaussian fit. In 
this momentum range, the tt", e~, K~ and p are clearly 
separated. The lower panel shows the dE/dx distribu- 
tion for the momentum range of 0.70 - 0.72 GeV/c. For 
the higher momentum bin, the tt^, e^ and K~ overlap; 
however, the p is still clearly separated. Even though 
it is still possible to extract kaons with 90% purity in 
this momentum bin, we only used the kaons up to 0.70 
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FIG. 2: Differential elliptic flow for pions, kaons and protons 
+ antiprotons for minimum-bias events. The solid lines show 
the fit with the modified Blast Wave model, and the dotted 
lines are the fit with the unmodified model. 



GeV/c. 

The flow analysis method [El| involves the calculation 
of the event plane angle, which is an experimental esti- 
mator of the real reaction-plane angle. For the analysis 
presented in this Letter, each particle was correlated with 
the event plane from all the other particles (for the other 
methods, see p2i). The differential elliptic flow, V2, de- 
pends on mass, rapidity (y) and pt- In Fig. H, V2{pt) 
is shown for pions, kaons, and protons -I- antiprotons 
for minimum-bias collisions 123] , integrated over rapidity 
and centrality by taking the multiplicity-weighted aver- 
age. The uncertainties shown are statistical only. Using 
the same procedure to estimate the systematic uncertain- 
ties as in Ref . [^ for V2 integrated over pt , the systematic 
uncertainty for minimum-bias data is 13%. We have ver- 
ified that the positive and negative identified particles 
used in this analysis have the same V2{pt) within statis- 
tical uncertainties. For V2{pt), the pions were integrated 
over \y\ < 1.0, the kaons over |y| < 0.8, and the protons 
+ antiprotons over \y\ < 0.5 (the rapidity ranges chosen 
correspond to approximately consistent jr/j coverage for 
all of the selected particles) . Mathematically the V2 value 
at Pt = 0, as well as its first derivative, must be zero. As 
a function of pi, the pions exhibit an almost linear depen- 
dence of V2 , whereas the protons -I- antiprotons exhibit a 
more quadratic behavior, clearly different from the pions. 
As expected in a hydrodynamic picture, the kaons lie be- 
tween the pions and the protons + antiprotons. The 90% 
purity of the protons + antiprotons in the 0.8 < pt < 0.9 
range leads to a maximum 4-5% systematic error on the 
V2 value in this bin. The observed behavior may be the 
result of the interplay between the mean expansion ve- 
locity, the elliptic component of the expansion velocity, 
and the thermal velocity of the particles. A similar effect 
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FIG. 3: Upper panel: Differential elliptic flow for pions in 
three different centrality bins. Lower panel: The same for 
protons -I- antiprotons. The dotted lines show the predictions 
from a full hydrodynamic model calculation [[l5|. The uncer- 
tainties shown are statistical only. 



was discussed for the case of directed flow in [g^ . 

The differential elliptic flow, W2(pt), is plotted for pions 
for three different centrality selections in the upper panel 
of Fig. ^, and for protons in the lower panel. The open 
triangles represent the most central 11% of the measured 
cross section. The open squares correspond to 11 - 45% 
of the measured cross section, and the open circles corre- 
spond to 45 - 85%. The uncertainties on the points are 
statistical only. The systematic uncertainty is smallest 
for the centrality region with the best reaction plane res- 
olution and is estimated to be 20% for the most central 
bin, 8% for the mid-central bin, and 22% for the most 
peripheral bin. At a given pt , the more peripheral colli- 
sions have the largest value oiv2{pt), and V2{pt) decreases 
for more central collisions. For all three of the centrality 
ranges, in the measured pt range, the pt dependence of 
V2 for pions is approximately linear. 

We have fitted the data with a simple hydrodynamic- 
motivated model. This model is a generalization of the 
Blast Wave model from ||l5|, 25 assuming the flow field 



is perpendicular to the freeze-out hyper-surface. 



V2{Pt) 
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r2TT 
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(1) 



where /q, /2, and Ki are modified Bessel functions. 



and where at((/)&) = (pt/T/) sinh(p(0&)), and /3t((^6) = 
{mt/Tf)cosh{p{(f>i,)). The basic assumptions of this 
model are boost-invariant longitudinal expansion and 
freeze-out at constant temperature Tf on a thin shell, 
which expands with a transverse rapidity exhibiting a 
second harmonic azimuthal modulation given by p{4'b) — 
Po + Pacos(20h). In this equation, 0(, is the azimuthal 
angle (measured with respect to the reaction plane) of 
the boost of the source element on the freeze-out hyper- 
surface |15[ , and po Etnd pa are the mean transverse ex- 
pansion rapidity {vq = tanh(po)) and the amplitude of 
its azimuthal variation, respectively. With S2 = 0, our 
equation reduces to Eq.2 from [g5|. In Fig. g, the fit to 
the minimum-bias data with S2 = is shown as the dot- 
ted lines. The poor fit shows that the data can not be 
described under the assumption of a spatially isotropic 
freeze-out hyper-surface in the transverse plane. This 
led us to generalize Eq.2 from ||l5[ to the case of a spa- 
tially anisotropic freeze-out hyper-surface, introducing 
one extra parameter, S2, describing the variation in the 
azimuthal density of the source elements, oc 2s2 cos(20fc). 
This additional parameter leads to a good description of 
the data, shown as the solid lines in Fig. g. A positive 
value of the S2 parameter would mean that there are more 
source elements moving in the direction of the reaction 
plane. The S2 parameter does not distinguish between 
an in-plane or out-of-plane extended source, because the 
direction of the flow flcld is not an observable. However, 
azimuthally sensitive HBT measurements will be able to 
address this. In this model (Eq. |l|), elliptic flow as a func- 
tion of particle transverse momentum, V2{pt), starts from 
zero and rises approximately quadratically with pt until 
the particle becomes relativistic, and then f2(pt) con- 
tinues to rise almost linearly. For the heavier particles, 
the linear rise is delayed. These mass-dependent effects 
are larger for lower temperatures (Tf) and larger trans- 
verse rapidities (po)- The linear rise increases with both 
the amplitude of the azimuthal variation of the trans- 
verse expansion rapidity (pa) and the elliptic deforma- 
tion (S2), but only pa produces a mass dependence. This 
arises from the dependence of the S2 parameter effect 
on momentum not energy. This also results in a change 
in slope and a characteristic curvature in V2{pt), since 
the flow associated with S2 saturates quite early, when 
Pt/T 3> 1 |£^. The fact that the data can not be de- 
scribed with S2 = leads to the interpretation that the 
elliptic flow is not caused by an azimuthal velocity varia- 
tion alone, but by the combination of a velocity difference 
and a spatially anisotropic freeze-out hyper-surface. 

Table | lists the results of the two flts (without, and 
with S2) for minimum-bias data. The parameters are 
correlated and therefore have large uncertainties. It is 
expected that better constraints on the parameters can 
be obtained in the future, when the flts to the single in- 
clusive spectra for different mass particles are available. 
As the solid lines in Fig. || show, the data could be de- 



scribed by a reasonable set of parameters. For the dif- 
ferent centralities, the Tf and po were kept fixed to the 
minimum-bias values. Both the obtained pa and S2, not 
shown, decrease from peripheral to central collisions as 
expected from a simple geometrical picture. 



r/(MeV) 


po 


Pa 


S2 


135 ± 19 
101 ± 24 


0.58 ± 0.03 
0.61 ± 0.05 


0.09 ± 0.02 
0.04 ± 0.01 




0.04 ± 0.01 



TABLE I: The parameters from the Blast Wave fit for 
minimum-bias collisions. The first row lists the parameters 
for the unmodified fit (s2 = 0), the second row gives the re- 
sults from the fit function including the spatially anisotropic 
freeze-out hyper-surface. 

A full hydrodynamic calculation, shown in Ref. p5[ , 
describes the v^ipt) for pions, kaons and protons -t- an- 
tiprotons for minimum-bias collisions equally well. The 
agreement is consistent with early local thermal equilib- 
rium and the presence of an early stage pressure gradient. 
However, within the current statistical uncertainties in 
the data and different theoretical interpretations ||2^ , the 
measurement does not justify drawing inferences about 
the different EOS. Comparing the centrality dependence 
from our Fig. Swith the same hydrodynamic calculation, 
shows a deviation which is most pronounced for the most 
peripheral data. This might indicate that the amount of 
particle rescattering in peripheral collisions is insufficient 
to justify the local thermal approximation implicit in hy- 
drodynamical models. 

We have made the first measurement of identified par- 
ticle elliptic flow at RHIC. The measured elliptic flow as 
a function of pt and centrality differs signiflcantly for par- 
ticles of different masses. This mass dependence can be 
described with a simple hydrodynamic-motivated Blast 
Wave model. This model suggests that elliptic flow is 
generated by the combination of an azimuthal velocity 
variation and a spatially anisotropic freeze-out hyper- 
surface. The mass dependence of V2{pt) is also in close 
agreement with full hydrodynamic model calculations, 
suggesting that in Au-|-Au coUisions at ,/s^ = 130 GeV, 
a system is created which, for central and mid-peripheral 
collisions, is consistent with early local thermal equilib- 
rium followed by hydrodynamic expansion. 
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